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Abstract 

The radiating part of a circular acoustic source is determined, on the 

P^ ' basis of an exact analysis of the radiation properties of a source with 

^-^, angular dependence expjn^ and arbitrary radial dependence. It is found 

' O ' that the number of degrees of freedom in the radiated field is no greater 

than k — n, where k is the wavenumber. The radiating part of the source at 
low frequency is explicitly stated and used to analyze noise cancellation. 
The results are applied to the identification of sources in jet noise and an 



^ 
^ 



Y^ I explanation for the low order structure of jet noise fields is proposed. 

^! 1 Introduction 

■ ^^' Sound generation by spinning modes is a central problem in many applications. 

Devices such as propellers and fans obviously produce a rotating source sys- 
tem, but the termination of a circular duct also radiates like a spinning source 
and jets have a source system which can be decomposed into spinning modes. 

^^ I The problem of the relationship between a source distribution and its spinning 

\^ . acoustic field has thus attracted considerable attention in the literature. 

This paper examines one part of the problem, the relationship between the 

i/^ . radial structure of the source and the form of the c orresponding a coustic field. 

CD ' Earlier work, using exact and asymptotic analysis (jCarlevl bOlOol U)!). has fixed 
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^D . an upper limit on the information which is radiated into the field, but without 

considering the effect of the source distribution. In this paper, the theory is 
extended to include the relationship between the source and the field, allowing 
a discussion of the implications for a number of problems. 

A first area where the analysis is relevant is that of control of noise from 
C^ ' rotors. One approach to this problem is to fit an inverse model to the measured 

noise and then use this model to com pute a noise field wh ic h canc els the rotor 
noise at some point. In a recent study ([Gerard et a/.l . l2005al la. l2007l) it was found 



that a very low order acoustic model of a cooling fan was sufficiently accurate 
for control purposes. At first glance, it is not clear why a low order model should 
give a good match to the results from a finite rotor, beyond considerations of 
acoustical compactness. An analysis of the information content of the field, 
however, shows that the field is generated by a set of low order modes with 




Figure 1: Coordinate system for disc radiation calculations 

the higher orde r modes being cut-off and generating exponentially small noise 
( Carlevl . lioTo^ '). 

The second application of the approach of this paper is to inverse methods, 
in which a 'source' is determined from acoustic field measurements. There have 
been numerous applications of such techniques, but of particular interest here is 
that of jet noise. It is known that the noise field o f a turbulent jet is re presented 



by a much lower order model than is the flo w (iJordan et all 120071) . This is 
partly explained by axial interference effects ( Freundl . l200l[ for example) but 
relatively little attention has bee n paid to the influence of the radial structure 
of the source until quite recently ( Michell . 120091) . 



This paper presents an analysis of the general problem of radiation from 
a disc source and fixes limits on the proportion of the source which actually 
radiates into the acoustic field, with no approximation other than the standard 
acoustic assumption of linearity. In particular, the far-field assumption is not 
required, making the results applicable over a wide range of parameters. The 
implications of the results are then discussed with respect to rotor and jet noise. 



2 Spinning sound fields 

The problem is formulated as that of calculating the acoustic field radiated by 
a monopole source distributed over a circular disc. This disc may be viewed 
as the source proper, such as in the case of rotor noise, or as part of a dis- 
tributed three-dimensional source, as in jet noise. The system for the analysis 
is shown in Figure [1] with cylindrical coordinates {r,6,z). All lengths are non- 
dimensionalized on disc radius and the disc lies in the plane ^ = 0. The field 
from one azimuthal mode of the acoustic source, specified as s„(ri) expjn^i, 
has the form Pn(fc, r, z) expjn^, with P„ given by the Rayleigh integral: 



P„(fc,r, z) = 
R' = 
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Figure 2: Transformation to equivalent line source 

where k is wavenumber and subscript 1 indicates variables of integration. The 
field due to higher order sources, such as dipoles and quadrupoles, would be 
found by differentiation of ([1]). The analysis to be presented does not include 
such sources but the conclusions drawn should still be valid. 



2.1 Equivalent line source expansion 

The analysis of the nature of the sound field from an arbitrary disc source 
is based on a transformation of the disc to an exactly equivalent line source, 
an approach which has been used to study tr ansient radiation from pistons 
( Oberhettingeij . Il961t IPiercd . 119891 ) , rotor noise dChapmanl . Il993t ICarlevi Il999l) 
and source identification methods (jCarlevl 120091) . 

The transformation to a line source is shown in Figure [H which shows the 
new coordinate system (r2,02:2) centred on a sideline of constant radius r. 
Under this transformation: 
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for observer positions with r > 1, with the limits of integration given by: 
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Functions of the form of K{r, r2) have been analyzed in previous work ( Carlevl 
19991 ). I n this paper, it is sufficient to note that the source function can be 
written (|Carlevl[2009[ ): 



K{r,r2) = 5]^.,(r)C/,(.)(l-s2)l/^ 

9=0 



(5) 



where Uq{s) is a Chebyshev polynomial of the second kind, s = r2 — r and the 
coefficients Uq{r) are functions of r but not of z. Inserting (O into ([2]): 



oo „i ]kR' 

P„(fc,r,z) = J2uq{r)J ^^^Uq{s)ir + s)il-sy/'ds, (6) 
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R' = [{r + .s)^ + zY ■ (7) 

The radiation pr operties of the integrals of ([6]) have been examined in some 



detail elsewhere (|Carlevl l2010d la) , giving an exact result for the in-plane case 
z = 0: 

P„(fc,r,0)=^eJ^'^f;u,(r)f (9 + 1)^^2^. (g) 

q=0 

For large order q, the Bessel function Jq{k) decays exponentially for fc < g so that 
the line source modes with order q > k are 'cut-off' and generate exponentially 
small noise fields. Since the integrals have their maximum in the plane z = 0, 
([S]) says that the whole field is exponentially small and modes Ug{r) with q > k 
are cut-off everywhere. This gives an indication of how much of a given source 
distribution radiates into the acoustic field, near or far. There only remains to 
establish the relationship between the radial source s„(ri) and the line source 
coefficients Uq{r). 

2.2 Series expansion for spinning sound fields 



A recently derived series (jCarlevl . l20I0ci) for the field radiated by a ring source 



of radius ri can be used to find an expression for the sound radiated by a disc 
source with arbitrary radial variation: 



|-27r gj(fci?'+nei) 
Rn ^ / dai , 

Jo 4^^' 



;2n^ 



TT I ^^_^^„(2n + 4m+I)(2m-I)!! 



4 (rii?)i/2 ^^^ > (2n + 2m)!! 

^ ^ m— 

Hi'i2™+i/2 (fc^)C+2™ (cos </.) J„+2™+i/2 (fcri ) , (9) 



with W-^'ix) the Hankel function of the first kind of order i/, Ji^ the Bessel func- 
tion of the first kind and P^ the associated Legendre function. The observer po- 
sition is specified in spherical polar coordinates R— [r^ -I- z^]^/^, (\) — tan~^ r j z. 
Multiplication by the radial source term ris„(ri) and integration gives an 
expression for the field radiated by a general source of unit radius and azimuthal 
order n: 

Pnik,r,z) = f--^l ± i-ir ^^" + '^^\f^; - ^^^ „V„(COS0)5„,2„, 

ni—0 ^ '" 
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Setting z == (0 = 7r/2, R^r) 



P„(fc,r,0) = -)^- 2-(2n + 2m)!! ^"^^'^^^^ 

771 — 

where use has been made of the expression ( Gradshtevn fc Rvzhikl . ll98Cl 8.756.1): 
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2.3 Line source coefRcients 



The expressions for P„ from §[2T2]and S I2.1l are both exact and can be equated to 
derive a system of equations relating the coefficients Uq{r) to the radial source 
distribution s„(ri): 
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Under repeated differentiation, (J12p becomes a lower triangular system of linear 
equations which connects the coefficients Uq{r) and Sn+2m'- 
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where superscript (u) denotes the wth partial derivative with respect to fc, eval- 
uated at fc = 0. 

Using standard series (JGradshtevn fc Rvzhikl . llQSOl) , the products of special 
functions can be written: 
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where [t/2] is the largest integer less than or equal to t/2, and 
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with the polynomials Vn.t and Wn,t given by: 
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Given the power series, the derivatives at fc = are readily found 
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Setting V = 0, 1, . . . yields an infinite lower triangular system of equations for 
Uq{r): 

EU = B, (18) 

with \J = [uq ui . . .]^ and the elements of matrix E and vector B given by: 
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Given a radial source term s„(ri), (IT51) can be solved to find the coefficients 
Ug(r) of the equivalent line source modes. Since it is lower triangular, the first 
few values of Ug can be reliably estimated, although ill-conditioning prevents 
accurate solution for arbitrary large q. 



2.4 Radiated field 



From the relationship between the radial source term and the line source coef- 
ficients, some general pr operties of the acoustic field can be stated. The first 
result, already found by ICarlevI ()2010al la) is that, since the line source modes 
with q + I > k are cut off, the acoustic field has no more than k degrees of 
freedom, in the sense that the radiated field is given by a weighted sum of the 
fields due to no more than k elementary sources. From pS|) . this result can be 
extended. 

The first extension comes from the fact that B2V+1 = 0, for v' < n, on the 
right hand side of ([T5| . This means that Ug, g = 2i;' + 1, is uniquely defined by 
the lower order coefficients with q < 2v' . The result is that the acoustic field 
of azimuthal order n, whatever might be its radial structure, has no more than 
k — n degrees of freedom. 

A second extension comes from examination of ([T5)) The first few entries of 
the system of equations are: 
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resulting in the solution: 
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so that the ratio of uq and ui is constant, for arbitrary s„(ri). This means that 
low frequency sources of the same radius and azimuthal order generate fields 
which vary only by a scaling factor, since the higher order terms are cut off. 

Finally, if we attempt to isolate a source s„ (ri ) associated with a single line 
source mode, by setting Uq = 1 for some q, with all other Uq = 0, we find that 



the line modes must occur in pairs, since if U2v' 
the condition B2v'+i = 0. 



1, U2v'+i 7^ 0, being fixed by 



3 Results 

To illustrate the application of the result of the previous section, we present some 
results for the calculation of the line source coefficients and for the use of the 
method to modify the radiating part of a source. We also discuss qualitatively 
the implications of the results for studies of jet noise. 



3.1 Line source coefficient evaluation 

The first results are a comparison of coefficients Uq{r) computed using ([T8|) and 
those computed directly from exact closed-form expressions for K{r, r2 ) in the 
case when the radial source term is a monomial in radius s„ = r" ( Carlev . 



19991 ). Figure [3] compares the two sets of coefficients for a ~ 0,2,4, with the 
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Figure 3: Line source mode coefficients computed using the method of S I2.3 
(sohd fines) and directly from analytical formulae (symbols) for r — 5/4, s — r^, 
a = (circles), a = 2 (squares) and a = 4 (diamonds) for n = 2, 4, 16. 
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Figure 4: Cancellation effects for radial source terms: top row radiated field from 
original s„(ri) (solid) and modified source s„(ri) — f3s[^{ri) (dashed); bottom 
row source terms s„(ri) (solid) and /3s'^{ri) dashed. The left hand column 
covers the case n = 2 and the right hand n = 8. In both cases, fc = 1, r = 5/4. 



plots terminated at a value of q where the difference between the two sets of 
results becomes noticeable, q « 20. This gives an indication of the effect of the 
ill conditioning of (IT51) . For q < 20, the computed values of Uq are reliable. It is 
noteworthy that for small q, the coefficients are practically equal for all values 
of a so that for low frequency radiation, the radiated fields will be practically 
indistinguishable . 



3.2 Radial cancellation 

The analysis so far has identified that part of a source distribution which radi- 
ates, as a function of wavenumber k. From the results, it appears that only a 
small part of the source is responsible for the acoustic field, with most of the 
line source modes being cut off throughout the field, except at high frequency. 
Indeed, asymptotic analysis (jCarlevl bOlOaI N) shows that even the cut-on modes 
radiate efficiently only into a small part of the acoustic field, with the exception 
of those of low order. In any case, this offers a method for examination of the 
radiation properties of a source. Given a source term s„(ri) the approach is to 
impose a secondary source s^(ri) which generates the same set of line source 
modes up to some required order. In the simplest case, we match the first line 



source coefficient uq, which automaticahy matches ui. If uq is known for both 
SniTi) and for s'^{ri), then the source s^(ri) — /3s^(ri) will have uq — ui — 0, if 
(3 is taken as the ratio of uq for the primary and secondary sources s„ (ri ) and 

Figure S] shows the results of such a procedure using sources s„(ri) = 
Jn{ari), with a the first extremum of J„ (similar to a duct mode), and s^ri) = 
1. The wavenumber fc = 1 and values n = 2, 8 have been used. In the first case, 
n = 2, the noise reduction in the plane z = is quite large, about 20dB, but 
there is a slight increase around z = I. This is because, as seen in Figure [31 
the coefficient U2 is quite large and is multiplied by a Bessel function of order 3, 
which is not of high enough order for the exponential decay with k which cuts 
off the mode. 

The cut-off behaviour is seen more clearly in the n = 8 case, where the 
reduction at z = is 40dB. At larger z, the redu ction is much s maller, but this 
is because, as found from asymptotic analysis (jCarlevl l2010ai r6l). the field in 



this region only contains contributions from the remaining lower order modes, 
starting with q = 2. 

3.3 Degrees of freedom in jet noise fields 

The results of i) l2.4[ regarding the number of degrees of freedom in the acoustic 
field, can help explain some features of experiments on jet noi se. Despite the lack 



of consensus on what is meant by the 'source' of jet noise ( Jordan fc Gervaisl 



2008t ISuzukil . |2010|) , some progress has been made by assuming that the source 



of jet noise can be identified with some combination of flow quantities. An open 
question, however, is which part of the source term radiates, since it is clear 
that only a small fraction of the flow generates the acoustic fleld. Two recent 
sets of res ults, one exper imental, the other numerical, illustrate the issues. 
In one. iFreundl ( 20011 ) has used direct numerical simulation to compute the 



flow and noise of a Mach 0. 9 jet, val i dating the noise prediction against exper- 



iments and showing that a iLighthilll (|1952[ ) source term accurately reproduces 



the acoustic fleld. Spatial filtering of the source, using a wavenumber criterion 
to remove the non-radiating part, left "a set of modes capable of radiating to 
the far field", with the caveat that "additional cancellation may occur due to 
the radial structure of the source which is not accounted for in this analysis" . 
Indeed, the radial st ructure of jet noise sources has not received much attention 



until quite recently ( Michell . 120091) 



An experimental result of some interest is that of [Jordan et ali (|2007l ) who 
performed a modal decomposition of a jet fiow field and a proper orthogonal 
decomposition optimized for the resolution of the far field noise. They found 
that more than 350 modes were needed to capture half of the flow energy while 24 
modes sufficed for 90% of the far-fleld noise. As they note, passage to the far 
field acts as a filter passing only a low-dimensional representation of the fiow. 

From these observations, it is plausible that the relatively low order structure 
of jet noise can be explained by the results of this paper. In the notation of this 
paper, k = irStM , where St is Strouhal number based on jet diameter and M 
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is jet M ach number. For the ran ge of Strouhal number important for jet noise 



g^ < 2 (iM ichalke fc Fuchsl . I1975D . k < 2ttM. For the M = 0.9 jet studied by 
Jordan et al. (2007|), for example, this yields k < 5.7 and no more than about 



six line source modes radiate from the axisymmetric source modes at the highest 
frequency of interest. This estimate would nee d to be modifi ed to take account 
of axial interference as in the far field analysis of lMichell ( 20091 ) but does offer the 



possibility of establishing some reasonable limits on the detail to be expected 
from acoustic measurements on jets and the requirements for low order models 
used in noise control. 



4 Conclusions 

The radiation properties of disc sources of arbitrary radial variation have been 
analyzed to establish the part of the source which radiates into the acoustic field, 
without recourse to a far field approximation. Limits have been established on 
the number of degrees of freedom of the part of the source which radiates and 
the implications of these limits have been discussed for the problems of rotor 
noise and studies of source mechanisms in jets. Future work will consider the 
use of the findings of this paper to study the radiating portion of full jet source 
distributions, including axial interference effects. 
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